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Valley currents and non-local resistances of graphene nanostructures with broken inversion symmetry are 
considered theoretically in the linear response regime. Scattering state wave functions of electrons entering the 
nanostructure from the contacts represented by groups of ideal leads are calculated by solving the Lippmann- 
Schwinger equation and are projected onto the valley state subspaces to obtain the valley velocity helds and 
total valley currents in the nanostructures. In the tunneling regime when the Fermi energy is in the spectral gap 
around the Dirac point energy, inversion symmetry breaking is found to result in strong enhancement of the non¬ 
local 4 terminal Btittiker-Landauer resistance and in valley currents several times stronger than the conventional 
electric current. These strong valley currents are the direct result of the injection of electrons from a contact into 
the graphene in the tunneling regime. They are chiral and occur near contacts from which electrons are injected 
into the nanostructure whether or not a net electric current flows through the contact. It is also pointed out that 
enhanced non-local resistances in the linear response regime are not a signature of valley currents arising from 
the combined effect of the electric held and Berry curvature on the velocities of electrons. 

PACS numbers: 72.80.Vp, 73.63.-b, 73.63.Rt 


I. INTRODUCTION 

Graphene is a single atomic layer of carbon atoms arranged 
on a honeycomb lattice. Since the early work of Wallace^ it 
has been known that the electronic energy bands of graphene 
near the Fermi energy take the form of a degenerate pair 
of cones, also referred to as ‘valleys.’ The electronic den¬ 
sity of states of graphene vanishes at the Dirac point energy 
which coincides with the Fermi energy in pristine graphene. 
More recently it was recognized that the electronic structure 
of graphene has topological properties that relate to the Berry 
phase and Berry curvature.^ Specifically, graphene has a non¬ 
zero Berry phase associated with closed paths in reciprocal 
space that enclose a valley apex.^ If the inversion symmetry of 
the graphene lattice is broken, for example, by a staggered po¬ 
tential at the two atoms of the unit cell, then the Berry curva¬ 
ture rik becomes non-vanishing near the apex of each valley."^ 
In semiclassical theories of electron transport in graphene, the 
electron velocity Vk is related to the Berry curvature by^’^“^ 

Vk = l^+kxf2k (1) 

h OK 

where Ck is the energy of a Bloch state with wave vector k 
and, in the absence of magnetic fields, 

h\i = qeEi ( 2 ) 

where qe is the electron charge and E is the electric field. 
Since Tik points in opposite directions in the two valleys, 
Eqs. 1 and 2 imply that, in the presence of an electric field 
E, electrons in the two valleys of graphene with broken inver¬ 
sion symmetry will have differing velocities. This difference 
in velocity might, in principle, be used to separate electrons 


belonging to the different valleys spatially, and thus be em¬ 
ployed in future valleytronic devices in which the valley in¬ 
dex of electrons plays a role somewhat analogous to that of 
the electron spin quantum number in spintronic devices. 

Recently, Gorbachev et al.^ carried out non-local 4-terminal 
resistance measurements on graphene placed on a hexagonal 
boron nitride (hBN) substrate, varying the graphene Fermi 
level by the application of a gate voltage. They found a strik¬ 
ing enhancement of the measured non-local resistance when 
the Fermi level passed through Dirac points for samples with 
the crystallographic axes of the graphene aligned with those 
of the hBN. They interpreted the enhanced non-local resis¬ 
tance as a signature of valley currents in their samples based 
on the following intuitive picture:^ The aligned hBN substrate 
breaks the inversion symmetry of graphene. This results 
in non-zero Berry curvatures for states close in energy to Dirac 
points. Then, according to Eqs. 1 and 2, an electric field E 
that drives electric current through the sample would result in 
differing electron velocities in the two graphene valleys. This 
would imply a net valley current transverse to the electric cur¬ 
rent when the Fermi level is close to a Dirac point energy.^ It 
was further argued^ that this valley current fiows into the re¬ 
gion of the sample between the voltage probes and induces an 
electric field there, resulting in a potential difference between 
the voltage probes and in an enhancement of the non-local re¬ 
sistance, consistent with the experiment.^ 

However, valley currents were not observed directly in the 
experiment of Gorbachev et al.^ Also, the interaction between 
the hBN and graphene is expected theoretically^^”^^ to open 
an energy gap in the electronic density of states of graphene 
around the Dirac point energy, where the Fermi level is lo¬ 
cated when the observed^ enhanced non-local resistance has 
its maximum value. The presence of this gap is supported by 
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the observation of activated transport by Gorbachev et al? in 
some of their samples. Within the gap the transport mecha¬ 
nism is quantum tunneling, a phenomenon that has no clas¬ 
sical analog. Therefore the applicability of the semiclassical 
Eqs. 1 and 2 to this regime is unclear and a fully quantum 
mechanical theory of transport is necessary. 

A rigorous, fully quantum mechanical approach to calcula¬ 
tions of multi-terminal transport coefficients of nanostructures 
is provided by Biittiker-Landauer theory. However, the con¬ 
ceptual framework of Biittiker-Landauer theory differs from 
that of the topological arguments^ that have been outlined 
above. For example, within Biittiker-Landauer theory, in the 
linear response regime (i.e., in the limit where the applied bias 
voltages and currents approach zero), the four-terminal resis¬ 
tances depend only on the set of quantum electron transmis¬ 
sion probabilities {Ti^j} between all pairs of contacts {i, j} 
and the refiection probabilities {Ri^i} at contacts {i} evalu¬ 
ated at the Fermi energy (or, at finite temperatures, at energies 
near the Fermi energy). In this limit, the electric field E that 
appears in Eq. 2 goes to zero, and therefore in this limit it 
has no effect on {Ti^j} and {Ri^i}. Therefore, according to 
Biittiker-Landauer theory, the topological mechanism of val¬ 
ley currents that is embodied in the term k x fik in Eq. 1 
has no effect on four-terminal resistances (or more generally 
on two-, three- or other multi-terminal resistances, including 
non-local resistances) in the linear response regime. Thus it 
follows from Biittiker-Landauer theory that non-local resis¬ 
tance measurements in the linear response regime can not pro¬ 
vide experimental evidence of topological effects arising from 
the electric field in Eq. 2. 

In other words, Biittiker-Landauer theory shows that non¬ 
local (and other) resistances measured in the linear response 
regime do not depend on whether electrons travel through the 
sample under the influence of an electric field due to applied 
bias voltages or simply scatter through the sample freely in the 
energy window between the highest and lowest contact elec¬ 
trochemical potentials without being subjected to any driving 
electric field. This means that the effects of the topological 
term k x fik in Eq. 1 that arise from the driving electric field 
cannot be detected by local or non-local resistance measure¬ 
ments in the linear response regime. 

It is therefore of interest to explore theoretically the multi¬ 
terminal resistances and valley currents in fully quantum me¬ 
chanical models of graphene devices with broken inversion 
symmetry from the perspective of Biittiker-Landauer theory. 
This is done in the present paper for transport in graphene 
nanostructures subjected to staggered potentials in the lin¬ 
ear response regime. It will be shown here that the applica¬ 
tion of a staggered potential results in strong enhancement of 
the non-local four-terminal resistance when the Fermi level 
is close to the Dirac point energy, in qualitative agreement 
with the experimental findings of Gorbachev et al.^ Despite 
the studied nanostructures having atomically abrupt bound¬ 
aries where electrons scatter strongly and crystal momentum 
is not conserved, valley currents, up to several times larger 
than the conventional electric current, will be shown to ap¬ 
pear in response to electrochemical potential differences be¬ 
tween the electrodes when the Fermi level is near the Dirac 


point of graphene nanostructures with broken inversion sym¬ 
metry. These large valley currents are not generated by the 
topological mechanism embodied in Eqs. 1 and 2 since the 
present calculations are in the linear response limit where the 
driving electric field E in Eq. 2 tends to zero and therefore 
it does not appear in the Hamiltonian of the system in these 
Biittiker-Landauer theory-based calculations. Because these 
strong valley currents occur in a gap in the energy spectrum 
of the nanostructure they require electron tunneling and con¬ 
sequently their strength decays rapidly as the distance from 
a contact that injects electrons into the graphene nanostruc¬ 
ture increases. These valley currents will be seen to be chiral 
and to travel along the edge of the graphene that is in con¬ 
tact with an electrode that injects electrons into the graphene. 
If electrons are injected into the graphene nanostructure with 
broken inversion symmetry from a scanning tunneling micro¬ 
scope (STM) tip the valley currents will be shown to form a 
vortex circulating around the location at which the electron in¬ 
jection occurs. At Fermi energies further from the Dirac point 
and outside of the gap in the density of states of the nanostruc¬ 
ture, valley currents are also induced by bias voltages applied 
to the nanostructure. However in this regime they are found 
to be weaker, not to require inversion symmetry breaking, and 
to extend into regions of the nanostructure that are not close 
to the contacts. 

The remainder of this paper is organized as follows. In Sec¬ 
tion II The model of graphene nanostructures with broken in¬ 
version symmetry coupled to current and voltage contacts that 
is studied in this work is presented. Biittiker-Landauer theory 
and the Lippmann-Schwinger equation and how they apply 
to this model are outlined. Valley currents, valley velocity 
fields and non-local resistances are defined and the methodol¬ 
ogy used to calculate them is described. The numerical results 
obtained from these calculations are presented in Section III. 
The significance of the present findings is discussed in Section 
IV. 


II. MODEL AND FORMALISM 

In this paper the graphene nanostructures will be described 
by the standard nearest neighbor tight-binding Hamiltonian on 
a honeycomb lattice, 

~ ^ ^ ^ ^ ^nm (^n^m h.C.^ , (3) 

n {n,m) 

where is the on-site energy, tnm = t = 2.7 eV defines 
the matrix element between pz orbitals on nearest-neighbor 
atoms and the spin index is suppressed. This Hamiltonian with 
Cn = 0 is known to describe the tt band dispersion of graphene 
well at low energies, and has been used in numerous studies 
of electron transport in graphene nanostructures.^^ In order 
to introduce inversion symmetry breaking into the model, the 
simple choice = ±A is made so that is positive on one 
atom of the graphene unit cell and negative on the other. The 
amplitude of the symmetry breaking energy is chosen to be 
A = 0.0602 eV consistent with estimates for graphene on 
hBN reported in Ref. 11 . 
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Figure 1. (Color online) Graphene nanostructure with armchair 
edges. The size of the nanostructure in the y-direction is 9.838 
nm. Electron stream is injected through current contact 1 and exits 
through current contact 2. There is no net electric current entering or 
leaving through the voltage contacts 1 and 2; the potential difference 
between them is measured. Wavy lines represent ideal semi-infinite 
ID leads connecting graphene C atoms to electron reservoirs. Upper 
right inset: The graphene sublattices A (B) = open (filled) circles. 
Lower right inset: Hexagonal (solid) and rhombic (dotted) Brillouin 
zones of graphene. K and K' are the two Dirac points. 

This idealized model has been chosen since the purpose 
of this paper is to investigate the fundamental effects of in¬ 
version symmetry breaking in its simplest form on multi¬ 
terminal transport coefficients and valley currents within the 
fully quantum mechanical Biittiker-Landauer framework. It 
should be noted that the graphene on hBN system is more in¬ 
volved since the lattice parameters of hBN and graphene have 
a 2% mismatch and local configurations with N atoms under 
the centers of graphene hexagons and B atoms under C atoms 
have the lowest energy. However, these complications will 
not be considered here. 

In Biittiker-Landauer theory, at zero temperature and in 
the linear response regime, the currents li flowing towards the 
nanostructure in contacts i are related to the electrochemical 
potentials of the contacts by 

Ii = j{Niin-Y,Ti,jH) (4) 

3 

where Ni is the number of electronic modes incident on the 
nanostructure from contact i and Ti^j is the multichannel 
electron transmission probability from contact j to contact i. 
Ti^i = Ri^i the multichannel electron reflection probability 
from the nanostructure in contact i. 

The coefficients Ti^j that enter the Biittiker-Landauer the¬ 
ory are calculated in this paper as in many previous theoretical 
studies of quantum transport in nanostructureswith 
semiconducting, molecular, metallic, 

magnetic,^^’^^’^^’^^ and carbon-based^^ constituents. Each 
contact is represented by a set of ideal semi-inflnite one¬ 
dimensional (ID) tight-binding leads (the wavy lines in Fig. 
1) with one orbital per site and nearest neighbor hopping. 


One such ideal lead is attached to each peripheral site of the 
graphene nanostructure that is adjacent to a contact. The 
Hamiltonian of lead n is 

H-Ln = E] - E] * 

r {r,s) 

where t is the same as in Hq^ (Eq. 3). The site energy in 
Eq. 5 is the same as that of the site of the graphene nanostruc¬ 
ture to which the lead is connected. The coupling Hamiltonian 
between lead n and the edge site of the graphene nanostruc¬ 
ture is 

Wn = + h.C.), ( 6 ) 

The quantum transmission amplitudes for an electron to scat¬ 
ter at energy E via the nanostructure from one ID ideal lead to 
another are found by solving the Lippmann-Schwinger equa¬ 
tion 

|^™) = |C) + Go(S)TV|^™), (7) 

where |0o^) is an electron eigenstate of the m}^ ideal semi- 
inflnite lead that is decoupled from the graphene nanostruc¬ 
ture, Go{E) is the Green’s function of the decoupled sys¬ 
tem of the ideal leads and the graphene nanostructure, and 
W = Wn is the coupling between the graphene nanos¬ 
tructure and the ideal leads. is the scattering eigenstate 
of the complete coupled system associated with the incident 
electron state |0o^). Then 

Tij{E) =J2\tiimiE)\^vi/vi (8) 

n,m 

where tl^^{E) is the quantum transmission amplitude [ob¬ 
tained from the scattering state deflned by Eq. (7)] for 
an electron at energy E to scatter via the graphene nanostruc¬ 
ture from ideal ID lead m of contact j to ideal ID lead n of 

contact i. The sum is over ideal leads n (m) in contact i (j). 

i(i) 

^n(m) electron velocity in ideal ID lead n (m) of contact 

i (j) at energy E\vl^ = where e are the energy eigenval¬ 
ues of the Hamiltonian i^Ln [Eq. 5] of an inflnite ideal ID 
tight binding chain. 

Having evaluated ^ at the Fermi energy in this way, the 
Biittiker equations 4 are solved in the linear response regime 
to And the non-local 4-terminal resistance 

i?NL = AU/J (9) 

where I is the current passing through the current contacts and 
AU = A/i/ge is the potential difference between the voltage 
contacts; see the contacts in Fig. 1. 

The valley currents induced in the nanostructure in response 
to to the electrochemical potential differences between the 
various contacts in the linear response regime are estimated 
as follows: The scattering state of electrons injected 

into the nanostructure from ideal ID lead m is calculated at 
the Fermi energy for every lead m by solving the Lippmann- 
Schwinger equation (7). Then the scattering state is pro¬ 
jected onto the subspaces of Bloch states of graphene that be¬ 
long to the K and K' valleys. This yields the projected valley 
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states, and respectively. For this purpose a Bloch 
state is assigned to the K (K') valley if its wave vector lies 
within the upper (lower) half of the rhombic Brillouin zone 
defined by the dotted boundary in the lower right inset of Fig. 
1 . 

The ^-component of the velocity operator for electrons 
within the graphene nanostructure is 

v^ = —[^,Hgn] ( 10 ) 

where ^ ^p^lap and is the ^—coordinate of atomic 

site p of the graphene nanostructure. The expectation value of 

in the graphene nanostructure in the state is then 

^ (n) 

^ k,l 

where k and I are nearest neighbor sites of the graphene 
nanostructure, = xi or yi, = {Zi\ 2 p^) and \Zi) is the 
2pz orbital of the carbon atom at site 1. 

For a graphene nanostructure with multiple contacts i each 
at its own electrochemical potential pi the electron transport 
through the device is governed by a weighted average over the 
velocities associated with the scattering states injected by the 
various ideal leads rrii that make up all of the contacts i. The 
relevant weighted average will be defined here as 

where Api = Pi—pmm and /imin is the electrochemical poten¬ 
tial of the contact with the lowest electrochemical potential. 

The weighted valley velocities and for electrons 
in valleys K and K' are defined similarly by replacing 
and in Eqs. (11) and (12) by their projections \'ip^) 

and and and onto the valleys K and K', 

respectively. The weighted valley velocity is then defined as 

V™' = (vk{ - VK'i) (13) 

Equation (11) expresses the expectation value of the elec¬ 
tron velocity as a sum of terms evaluated at pairs (/c, /) of 
nearest neighbor atoms of the graphene nanostructure. This 
suggests that each such term be interpreted as the value of 
the velocity field for the scattering state at the loca¬ 

tion of each nearest neighbor pair. Assigning this value of 
the velocity field arbitrarily to the mid point {x^y) = 
{{xk + xi)/2, {yk + yi)/2) of the atomic pair yields 

-6)(or (i4) 

The corresponding velocity field weighted by the contribu¬ 
tions of the states contributing to transport through the nanos¬ 
tructure is then, as in Eq. (12), given by 

vl{x,y) = A/Ui (15) 

m,i m,i 
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Figure 2. (Color online) Calculated linear response properties vs. 
Fermi energy of the structure in Fig. 1 at zero temperature. Cur¬ 
rent I flows through the current contacts with no net current through 
either voltage contact. Results for symmetry breaking parameter 
A = 0.0602 (0.0) eV are orange (black), (a) Nonlocal resistance 
i^NL [Eq. (9)]. (b) Normalized valley velocities and vl^^jvy 

[Eq. (12), (13)] are solid and dashed lines, respectively. 

Obtaining the weighted velocity fields and for the 
valleys K and K' similarly by replacing by the projec¬ 
tions and the valley velocity field is defined as 

v'f^{x,y) = v^^{x,y) - v^,^{x,y) (16) 


III. RESULTS 

The results of Biittiker-Landauer calculations of the non¬ 
local four-terminal resistance i?NL defined by Eq. (9) for the 
structure in Fig. 1 at zero temperature in the linear response 
regime are shown in Fig. 2(a) as a function of the Fermi 
energy The Fermi level crosses the Dirac point energy 
at E’f = 0. Near the Dirac point, i?NL foi* the model with 
A = 0.0602eV ( orange line) exceeds i?NL for A = 0 (black 
line) by a factor of ^ 2.5. Thus within Biittiker-Landauer the¬ 
ory, the breaking of inversion symmetry of the graphene unit 
cell results in strong enhancement of the non-local resistance 
near the Dirac point. However, as can also be seen in Fig. 2(a), 
well away from the Dirac point energy the inversion symmetry 
breaking has little effect on the non-local resistance. 

Since the net electron fiow in the structure in Fig. 1 is 
from current contact 1 to current contact 2 (i.e., in the y- 
direction), the weighted velocity vector v (Eq. 12) is expected 
to point in the ^-direction. Accordingly, Vx in Eq. 12 is 
found to be zero within numerical error in the present compu¬ 
tations. The computed normalized, weighted valley velocities 
jvy and jvy are shown in Eig. 2(b). They are found 
to be non-zero except at isolated values of the Eermi energy 
both in the presence and absence of symmetry breaking. Note 
that /vy = /I where P^^ is the ^-component of the 
valley current and / is the total electric current through the 
nanostructure. The most striking feature of Eig. 2(b) is the 
strong peak near the Dirac point (E’f = 0) of v^^^jvy = Px^/1 
for broken inversion symmetry, (the solid orange curve). At 
its maximum the valley current in the x-direction exceeds the 
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current contact 1 voltage contact 1 


Figure 3. (Color online) Valley velocity field [Eq. (16)] in the lower part of the structure in Fig. 1 for electron flow from current contact 1 to 
current contact 2 (shown in Fig. 1). For electron fiow in this direction, the valley velocity field is much stronger in the vicinity of the contacts 
shown here than elsewhere in the graphene nanostructure. A = 0.0602eV. Ey — 0. 


total conventional electric current through the nanostructure 
by a factor of more than 4. This peak of the valley current 
is in the gap (of width 0.264 eV) in the energy spectrum of 
the broken symmetry nanostructure around the Dirac point. 
By contrast, the x-component of the valley current in the ab¬ 
sence of symmetry breaking (solid black curve) and the y- 
component of the valley current with (orange dashed curve) 
and without (black dashed curve) symmetry breaking all van¬ 
ish at the Dirac point and are relatively weak elsewhere. As 
can also be seen in Fig. 2(b), for well away from the Dirac 
point energy, both and Vy^^jvy are insensitive to the 

breaking of the inversion symmetry of the graphene. 

The full widths of the central peaks of both the non-local 
resistance and the normalized valley velocity v^J^^jvy for the 
system with broken inversion symmetry [the orange curves 
in Figs. 2(a) and (b), respectively] are close in size to the 
0.264eV gap in the energy spectrum of the broken symmetry 
nanostructure. Because the value of the symmetry breaking 
parameter is relatively small (A = 0.0602eV), the size of the 
spectral gap is determined mainly by the quantum confine¬ 
ment of the electrons in the graphene nanostructure and the 
armchair character of the nanostructure’s edges. Thus for the 
same nanostructure but with the symmetry breaking turned off 
(A = 0) the width of the energy gap has a similar value, 0.235 
eV. For this reason the main nonlocal resistance peak in Fig. 
2(a) has almost the same width for A = 0 (the black curve) as 
for A = 0.0602eV (the orange curve). 

The valley velocity field y) is shown in Fig. 3 for 

the lower part of the nanostructure in Fig. 1. Inversion sym¬ 
metry is broken and the Fermi level is at the Dirac point. The 
valley velocity is large near current contact 1 and voltage con¬ 
tact 1. Its magnitude initially increases but then decreases 
rapidly with increasing distance from the contacts. The val¬ 
ley velocity is clearly chiral, pointing mainly from left to right 
near the graphene boundary shown (its overall direction re¬ 
verses if the sign of A, the symmetry breaking parameter, is 
changed) but it does not extend along the boundary much be¬ 
yond where a contact ends. 

In Fig. 3 the electron flow enters the graphene nanostruc¬ 
ture through current contact 1 and exits through current con¬ 
tact 2 that is located outside of the region shown in Fig. 3; 
see Fig. 1 for its location. If the direction of the electron flow 
through the graphene nanostructure is reversed, so that elec- 




Figure 4. (Color online) (a) Valley velocity field [Eq. (16)] and 
(b) velocity field [Eq. (15)] for electrons injected into a graphene 
nanostructure with broken inversion symmetry via a single carbon 
atom (orange). Only a small part of the graphene nanostructure is 
shown. A = 0.0602eV. E'f = 0. 


trons flow instead from current contact 2 to current contact 1 
(in Fig. 1), then the valley velocity fleld becomes strongest 
near current contact 2 and voltage contact 2, i.e., near the op¬ 
posite edge of the sample to that where the valley velocity fleld 
is strongest in Fig. 3. For electron flow from current contact 
2 to current contact 1 the direction of the valley velocity fleld 
is from right to left, i.e., its direction is opposite to that in Fig. 
3, consistent with the chiral character of the valley current. 

The chiral nature of the valley current is further clarifled in 
Fig. 4. Fig. 4 (a) shows the strongest part of the valley ve¬ 
locity fleld associated with electron injection into a graphene 
nanostructure with broken inversion symmetry via a single in¬ 
terior carbon atom (colored orange) of the nanostructure, as 
in an idealized STM setup. The valley velocity fleld forms a 
vortex circulating clockwise (counter-clockwise if the sign of 
A is changed) around the injection point, whereas the electron 
flux travels outwards overall from the injection point as shown 
in Fig. 4 (b). 

In the vicinity of a contact the valley current is due almost 
entirely to electrons injected into the nanostructure from that 
contact. This is true even for the voltage contact (Fig. 3) 
through which no net electric current flows since the zero net 
current is due to equal fluxes of electrons entering and leav¬ 
ing the contact. The valley current of electrons leaving the 
contact is much larger than that of those entering the contact. 
Thus a contact through which no net electric current flows can 
be used to create a valley current into a graphene nanostruc¬ 
ture with broken inversion symmetry. I.e., it can in principle 
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generate a pure valley current. 

The results presented above have been for a graphene 
nanostructure with only armchair edges. Similar calcula¬ 
tions for a rectangular structure of similar size (dimensions 
9.656nm x 9.838nm) and similar contacts but with two zigzag 
and two armchair edges were carried out and yielded quali¬ 
tatively similar results but the effects of inversion symmetry 
breaking were found to be much stronger in this case: The 
non-local resistance for A = 0.0602eV was found to exceed 
that for A = 0 by more that a factor of 100 for E’p near 
the Dirac point energy. Also, at its maximum the valley cur¬ 
rent in the x-direction was found to exceed the total conven¬ 
tional electric current through the nanostructure by a factor of 
more than 18. These large numbers are attributable to the flat 
electronic dispersion at graphene zigzag edges^^ in tight bind¬ 
ing models described by the non-interacting electron Hamilto¬ 
nian, Eq. 3. However, theoretical studies have suggested that 
electron-electron interactions may give rise to magnetism at 
zigzag edges.The potential implications of this for valley 
currents and nonlocal resistances are beyond the scope of the 
present paper. 

Different approaches for realizing valley currents in 
graphene have also been proposed based on electric flelds act¬ 
ing on electrons in the presence of Berry curvature,"^ graphene 
point contacts with zigzag edges,electron scattering at 
the boundary between monolayer and bilayer graphene,^^ 
electron scattering at a line defect in graphene,illumina¬ 
tion of monolayer^^ or bilayer^^ graphene by circularly po¬ 
larized radiation, optical injection of a pure valley current 
in graphene,^^ and gate-induced valley Altering in bilayer 
graphene."^^ Strong valley current polarizations, comparable 
to those obtained with the present approach for the structure 
having both zigzag and armchair edges have been estimated 
for some of these approaches,^"^’^^’^^""^^ albeit for models of in- 
flnite two-dimensional graphene or inflnite graphene ribbons. 
Also, unlike in the present work, the spatial distribution of val¬ 
ley currents was not reported in the previous studies,^’^"^”^^’^^ 
and the effects of boundaries between the graphene and source 
and drain electrodes were not taken into account. 


IV. CONCLUSIONS 

The present work suggests an approach for creating valley 
currents in graphene with broken inversion symmetry that dif¬ 


fers fundamentally from previous proposals.As has 
been explained above, it follows from Biittiker-Landauer the¬ 
ory that in the linear response regime considered here, the 
transport properties of nanostructures are determined by elec¬ 
tron scattering states calculated in the limit where the driving 
electric field has been sent to zero. Consequently in the lin¬ 
ear response regime the acceleration of electrons by the driv¬ 
ing electric held has no effect on multi-terminal resistance 
coefficients or on jl, the ratio of the valley current 
and the conventional electric current / passing through the 
nanostructure. Thus the valley currents discussed here are not 
due to electron acceleration in an electric field in the pres¬ 
ence of Berry curvature but instead are a direct consequence 
of non-adiabatic injection of electrons from a contact into the 
graphene. They are strongest for graphene with broken in¬ 
version symmetry in the tunneling regime when the Fermi 
energy is in the spectral energy gap around the Dirac point. 
Consequently in this regime these valley currents are strongest 
close to the graphene/contact boundary. They are chiral and 
can be very strong close to the Dirac point, i.e., several times 
larger than the conventional electric current even after averag¬ 
ing over the entire graphene nanostructure. They can appear 
even at a voltage contact through which no net conventional 
electric current flows provided that electrons are being emitted 
(and absorbed) by that contact. They are predicted to be real¬ 
ized whenever electrons cross into the graphene at an abrupt 
boundary (which may be regular or rough on the atomic scale) 
between a contact and graphene with broken inversion sym¬ 
metry at energies in the spectral gap around the Dirac point. 
At Fermi energies well away from the Dirac point (i.e., out¬ 
side of the gap in the density of states of the nanostructure) 
valley currents can still be induced by bias voltages applied 
to the nanostructure but in this regime they are considerably 
weaker, are not sensitive to whether or not the inversion sym¬ 
metry of the graphene is broken and are not confined to re¬ 
gions of the nanostructure that are close to contacts. That val¬ 
ley currents can be induced in graphene nanostructures even 
in the absence of inversion symmetry breaking has been rec¬ 
ognized previously,and is a consequence of the fact that 
the Bloch state wave vector need not be a conserved quantity 
in nanostructures whose translational crystal symmetries are 
broken due to the presence of boundaries. 

This research was supported by NSERC, CIFAR, Westgrid, 
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